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Mass bounds on dark matter (DM) candidates are obtained for particles that decouple in or 
out of equilibrium while ultrarelativistic with arbitrary isotropic and homogeneous distribution 
functions. A coarse grained Liouville invariant primordial phase space density T> is introduced which 
depends solely on the distribution function at decoupling. The density T> is explicitly computed and 
combined with recent photometric and kinematic data on dwarf spheroidal satellite galaxies in the 
Milky Way (dShps) and the observed DM density today yielding upper and lower bounds on the 
mass, primordial phase space densities and velocity dispersion of the DM candidates. Combining 
these constraints with recent results from TV-body simulations yield estimates for the mass of the 
DM particles in the range of a few keV. We establish in this way a direct connection between the 
microphysics of decoupling in or out of equilibrium and the constraints that the particles must fulfill 
to be suitable DM candidates. If chemical freeze out occurs before thermal decoupling, light bosonic 
particles can Bose-condense. We study such Bose-Einstein condensate (BEC) as a dark matter 
candidate. It is shown that depending on the relation between the critical (T c ) and decoupling (To) 
temperatures, a BEC light relic could act as CDM but the decoupling scale must be higher than the 
electroweak scale. The condensate hastens the onset of the non-relativistic regime and tightens the 
upper bound on the particle's mass. A non-equilibrium scenario which describes particle production 
and partial thermalization, sterile neutrinos produced out of equilibrium and other DM models is 
analyzed in detail and the respective bounds on mass, primordial phase space density and velocity 
dispersion are obtained. Thermal relics with m ~ fewkeV that decouple when ultrarelativistic and 
sterile neutrinos produced resonantly or non-resonantly lead to a primordial phase space density 
compatible with cored dShps and disfavor cusped satellites. Light Bose-condensed DM candidates 
yield phase space densities consistent with cores and if T c 3> Id also with cusps. Phase space 
density bounds on particles that decoupled non-relativistically combined with recent results from 
N-body simulations suggest a potential tension for WIMPs with m ~ 100 GeV,Td ~ lOMeV. 
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I. INTRODUCTION 

Although the existence of dark matter (DM) was inferred several decades ago its nature still remains elusive. 
Candidate dark matter particles are broadly characterized as cold, hot or warm depending on their velocity dispersions. 
The clustering properties of collisionless DM candidates in the linear regime depend on the free streaming length, 
which roughly corresponds to the Jeans length with the particle's velocity dispersion replacing the speed of sound in 
the gas. Cold DM (CDM) candidates feature a small free streaming length favoring a bottom-up hierarchical approach 
to structure formation, smaller structures form first and mergers lead to clustering on the larger scales. 

Among the CDM candidates are weakly interacting massive particles (WIMPs) with m ~ 10 — 10 2 GeV. Hot 
DM (HDM) candidates feature large free streaming lengths and favor top down structure formation, where larger 
structures form first and fragment. HDM particle candidates are deemed to have masses in the few eV range, and 
warm DM (WDM) candidates are intermediate with a typical mass range m ~ 1 — lOkeV. 

The concordance ACDM standard cosmological model emerging from CMB, large scale structure observations and 
simulations favors the hypothesis that DM is composed of primordial particles which are cold and collisionless 0. 
However, recent observations hint at possible discrepancies with the predictions of the ACDM concordance model: 
the satellite and cuspy halo problems. 



The satellite problem, stems from the fact that CDM favors the presence of substructure: much of the CDM is 
not smoothly distributed but is concentrated in small lumps, in particular in dwarf galaxies for which there is scant 
observational evidence so far. A low number of satellites have been observed in Milky- Way sized galaxies 0, 0, 
This substructure is a consequence of the CDM power spectrum which favors small scales becoming non-linear first, 
collapsing in the bottom-up hierarchical manner and surviving the mergers as dense clumps [1, Q . 

The cuspy halo problem arises from the result of large scale A-body simulations of CDM clustering which predict 
a monotonic increase of the density towards the center of the halos 0, H, H, E3|, for example the universal 
Navarro-Frenk- White (NFW) profile p(r) ~ r _1 (r + tq)~ 2 Q which describes accurately clusters of galaxies, but 
indicates a divergent cusp at the center of the halo. Recent observations seem to indicate central cores in dwarf 
galaxies [HI H2, EH, EH , leading to the 'cusps vs cores' controversy. 

A recent compilation of observations of dwarf spheroidal galaxies dSphs [l4j, which are considered to be prime 
candidates for DM subtructure [l5j . seem to favor a core with a smoother central density and a low mean mass 
density ~ 0.1M Q /pc 3 rather than a cusp [II]. The data cannot yet rule out cuspy density profiles which allow a 
maximum density < 60M Q /pc 3 and the interpretation and analysis of the observations is not yet conclusive [ITI.[l^|. 
These possible discrepancies have rekindled an interest in WDM particles, which feature a velocity dispersion larger 
than CDM particles, and consequently larger free-streaming lengths which smooth-out the inner cores and would be 
prime candidates to relieve the cuspy halo and satellite problems (l7| . 

A possible WDM candidate is a sterile neutrino [H, ES H3| with a mass in the keV range and produced via their 
mixing and oscillation with an active neutrino species either non- resonantly [l8| . or through MSW (Mikheiev-Smirnov- 
Wolfenstcin) resonances in the medium [l9| . Sterile neutrinos can decay into a photon and an active neutrino (more 
precisely the largest mass eigenstate decays into the lowest one and a photo n) 1 211) yielding the possibility of direct 
constraints on the mass and mixing angle from the diffuse X-ray background [22j . 

Observations of cosmological structure formation via the Lyman-a forest provide a complementary probe of primor- 
dial density fluctuations on small scales which yield an indirect constraint on the masses of WDM candidates. While 
constraints from the diffuse X-ray background yield an upper bound on the mass of a putative sterile neutrino in the 
range 3 — 8keV [22j], the latest Lyman-a analysis [23j yields lower bounds in the range 10 — 13keV in tension with 
the X-ray constraints. More recent constraints from Lyman-a yield a lower limit for the mass of a WDM candidate 
TnwDM ^ 1.2keV(2tr) for an early decoupled thermal relic and tuwdm ^ 5.6keV(2<r) for sterile neutrinos [2~jj . 
Strong upper limits on the mass and mixing angles of sterile neutrinos have been recently discussed [25(, however, 
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there are uncertainties as to whether WDM candidates can explain large cores in dSphs [26|. It has been recently 
argued [27j that if sterile neutrinos are produced non-resonantly [l8[ the combined X-ray and Lyman-a data suggest 
that these cannot be the only WDM component, with an upper limit for their fractional relic abundance < 0.7. Recent 
[HI constraints on a radiatively decaying DM particle from the EPIC spectra of (M31) by XMM-Newton confirms 
this result and places a stronger lower mass limit m < 4 keV. 

All these results suggest that DM could be a mixture of several components with sterile neutrinos as viable candi- 
dates. 

Motivation and goals: Although the ACDM paradigm describes large scale structure formation remarkably well, 
the possible small scale discrepancies mentioned above motivate us to study new constraints that different dark matter 
components must fulfill to be suitable candidates. Cosmological bounds on dark matter components primarily focused 
on standard model neutrinos [29|, [3(| , heavy relics that decoupled in local thermodynamic equilibrium (LTE) when 
non-relativistic [U [H, HU or thermal ultrarelativistic relics [HI, HE HI, I37ll3pl. More recently, cosmological precision 
data were used to constrain the (HDM) abundance of low mass particles |39l.l40l. l4ll [42j assuming these to be thermal 
relics. 

The main results of this article are: 

(a:) We consider particles that decouple in or out of LTE during the radiation dominated era with an arbitrary 
(but homogeneous and isotropic) distribution function. Particles which decouple being ultrarelativistic eventually 
become non-relativistic because of redshift of physical momentum. We establish a direct connection between the 
microphysics of decoupling in or out of LTE and the constraints that the particles must fulfill to be suitable DM 
candidates in terms of the distribution functions at decoupling. 

(b:) We introduce a primordial coarse grained phase space density 

n(t) 



V 



where n(t) is the number of particles per unit physical volume and (Pf) i s the average of the physical momentum 
with the distribution function of the decoupled particle. I? is a Liouville invariant after decoupling and only depends 
on the distribution functions at decoupling. In the non-relativistic regime T> is simply related to the phase densities 
considered in refs. [ll|, H3, US HH] and can only decrease by collisionless phase mixing or self-gravity dynamics [43| . 
In the non-relativistic regime we obtain 

v = ^hT~ (I'D 

3= m 4 v DM 

where a dm is the primordial one-dimensional velocity dispersion and pdm the dark matter density. Combining the 
result for the primordial phase space density T> determined by the mass and the distribution function of the decoupled 
particles, with the recent compilation of photometric and kinematic data on dSphs satellites in the Milky- Way [3| 
yields lower bounds on the DM particle mass m whereas upper bounds on the DM mass are obtained using the 
value of the observed dark matter density today. Therefore the combined analysis of observational data from (dSphs), 
N-body simulations and the present DM density allows us to establish both upper and lower bounds on the mass of 
the DM candidates. 

We thus provide a link between the microphysics of decoupling, the observational aspects of dark matter halos and 
the DM mass value. 

(c:) Recent A-body simulations Q indicate that the phase-space density decreases a factor ~ 10 2 during gravita- 
tional clustering. This result combined with eq. fll.ip and the observed values on dSphs satellites [HI yield 

2 8 

Scored ^ T ^eV , rn CUS p <~ — keV . 

9 4 3 4 

for the masses of thermal relics DM candidates, where 'cored' and 'cusp' refer to the type of profile used in the dShps 
description and 1 < g < 4 is the number of internal degrees of freedom of the DM particle. Wimps with masses 
~ 100 GeV decoupling in LTE at temperatures ~ 10 MeV lead to primordial phase space densities many orders of 
magnitude larger than those observed in (dSphs). The results of A-body simulations, which yield relaxation by 2 — 3 
orders of magnitude [44] suggest a potential tension for WIMPs as DM candidates. However, the A-body simulations 
in ref.[44j begin with initial conditions with values of the phase space density much lower than the primordial one. 
Hence it becomes an important question whether the enormous relaxation required from the primordial values to 
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those of observed in dSphs can be inferred from numerical studies with suitable (much larger) initial values of the 
phase space density. 

(d:) We study the possibility that the DM particle is a light Boson that undergoes Bose-Einstein Condensation 
(BEC) prior to decoupling while still ultrarelativistic. (This possibility was addressed in [35[). We analyze in detail 
the constraints on such BEC DM candidate from velocity dispersion and phase space arguments, and contrast the 
BEC DM properties to those of the hot or warm thermal relics. 

(e:) Non-equilibrium scenarios that describe various possible WDM candidates are studied in detail. These scenarios 
describe particle production j4|| and incomplete thermalization (46j, resonant [l8| and non-resonant [n| production 
of sterile neutrinos and a model recently proposed [2(| to describe cores in dSphs. 

Our analysis of the DM candidates is based on their masses, statistics and properties at decoupling (being it in 
LTE or not). We combine observations on dSphs [l4| and TV-body simulations [44j |. with theoretical analysis using 
the non- increasing property of the phase space density [HI, [3(| [H, |43[ • 

The results from the combined analysis of the primordial phase space densities, the observational data on dSphs 
[l4j and the TV-body simulations in ref.[44[ are the following: 

• (i): conventional thermal relics, and sterile neutrinos produced resonantly or non-resonantly with mass in the 
range m ~ fewkeV that decouple when ultrarelativistic lead to a primordial phase space density of the same 
order of magnitude as in cored dShps and disfavor cusped satellites for which the data [l4| yields a much larger 
phase space density. 

• (ii): CDM from wimps that decouple when non-relativistic with m > 100 GeV and kinetic decoupling at 
T d ~ 10 MeV [H| yield phase space densities at least eighteen to fifteen orders of magnitude [see eqs. (|4.30[) . 
(|4.31[) and (|4.37|) ] larger than the typical average in dSphs ,14j. Results from TV-body simulations, albeit 
with initial conditions with much smaller values of the phase space density, yield a dynamical relaxation by 
a factor 10 2 — 10 3 [44|. If these results are confirmed by simulations with larger initial values there may be 
a potential tension between the primordial phase space density for thermal relics in the form of WIMPs with 
m ~ 100 GeV, T d ~ 10 MeV and those observed in dShps. 

• (iii) : Light bosonic particles decoupled while ultrarelativistic and which form a BEC lead to phase space densities 
consistent with cores and also consistent with cusps if T c /Td > 10. However if these thermal relics satisfy the 
observational bounds, they must decouple when gd g~ J (Td/T c )s > 130, namely above the electroweak scale. 

Section II analyzes the generic dynamics of decoupled particles for any distribution function, with or without LTE 
at decoupling, and for different species of particles. In section III we consider light thermal relics which decoupled 
in LTE as DM components: fermions and bosons, including the possibility of a Bose-Einstein condensate. Section IV 
deals with coarse grained phase space densities which arc Liouville invariant and the new bounds obtained with them 
by using the observational dSphs data and recent results from TV-body simulations, bounds from velocity dispersion, 
and the generalized Gunn-Trcmaine bound. In Section V we study the case of particles that decoupled out of 
equilibrium and the consequences on the dark matter constraints. Section VI summarizes our conclusions. 



II. PRELIMINARIES: DYNAMICS OF DECOUPLED PARTICLES 

While the study of kinetics in the early Universe is available in the literature [32], [4?], [48[ , in this section we expand 
on the dynamics of decoupled particles emphasizing several aspects relevant to the analysis that follows in the next 
sections. 

Consider a spatially flat FRW cosmology with length element 

ds 2 = dt 2 - a 2 (t) dx 2 (2.1) 

the non-vanishing Christoffcl symbols are 

T% = da5 t3 , r;„ (2.2) 

The (contravariant) four momentum is defined as p^ — dx^/dX with A an affine parameter, so that g tiu p^p u = m , 
where m is the mass of the particle. This leads to the dispersion relation 

p°(t) = Vm 2 + a 2 (t) £2(i) . (2.3) 
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The geodesic equations are 



JJ=-2#(i)p (i)f(i) p=-2H(t)p[t), (2.4) 

where = f and we used d/rfA = p° d/dt. The solution of eq. (|2.4p is 

where p c is the time independent comoving momentum. The local observables, energy and momentum as measured 
by an observer at rest in the expanding cosmology are given by 

E(t) = e% p» , P}(i) = -<^ef^ (2.6) 

where e£ form a local orthonormal tetrad (vierbein) 

9nv £q = ?7q/3 = diag(l, -1, -1, -1) , 

and the sign in eq. (|2.6[) corresponds to a space-like component. For the FRW metric 



(2.7) 



and we find. 



E=p a , P f ( t ) = a (t) p{t) = ^ ■ (2-8) 

Pf is clearly the physical momentum, redshifting with the expansion. Combining the above with eq. (|2.3p yields the 
local dispersion relation 

E{t)=p\t) = v /m 2 +P / 2 W. (2.9) 

A frozen distribution describing a particle that has been decoupled from the plasma is constant along geodesies, 
therefore, taking the distribution to be a function of the physical momentum Pf and time, it obeys the Liouville 
equation or collisionlcss Boltzmann equation 

* «lg*-o. (2.0, 

Taking Pf as an independent variable this equation leads to the familiar form 
Obviously a solution of this equation is 

f[P r ,t} = f d [a(t)P f }=f d \p c ] , (2.12) 

where p c is the time independent comoving momentum. The physical phase space volume element is invariant, 
d 3 Xf d 3 xPf = d 3 x c d 3 p c , where /, c refer to physical and comoving volumes respectively. 
The scale factor is normalized so that 

, . 1 + za , 

a « = TTW) (2 ' 13) 

and Pf(td) = Pci where t& is the cosmic time at decoupling and z is the redshift. 
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If a particle of mass m has been in LTE but it decoupled from the plasma with decoupling temperature T4 its 
distribution function is 

f d ( Pc ) = _4 — > ( 2 - w ) 

e T <* ±1 

for fermions (+) or bosons (— ) respectively allowing for a chemical potential pd at decoupling. 
In what follows we consider general distributions as in eq. (|2.12p unless specifically stated. 
The kinetic energy momentum tensor associated with this frozen distribution is given by 

T » = 9 J ^ if fdiPc) ' (2 ' 15) 

where g is the number of internal degrees of freedom, typically 1 < g < 4. Taking the distribution function to be 
isotropic it follows that 

d 3 P f 



r flop 

TS = 9 J^P°MPo)=P (2.16) 



J 3 

where p is the energy density and V is the pressure. In summary, 



fdHt) pa ; v = 1 y w 7=fe /dK<) F/] • (2 - i8) 



The pressure can be written in a manner more familiar from kinetic theory as 



V 



d 3 P f \v f \ 2 
(2tt) 3 



+ P] f d [a(t) P f ] , (2.19) 



where Vf — Pf/E is the physical (group) velocity of the particles measured by an observer at rest in the expanding 
cosmology. 

To confirm covariant energy conservation recall that d 3 Pf = d 3 p c /a 3 (t); Pf — p c /a(t); dfd/dt = 0, furthermore 
from cq.([23) it follows that p° = -H(t) Pf/p°, leading to 



d 3 P f P[ 
(2tt) 3 p° 



-3 H(t) p - H(t) g j — L -± f d ( Pc ) , ( 2 .20) 



the first term results from the measure and the last term from p°; from the expression of the pressure eq. (|2.17p the 
covariant conservation equation 

p + 3H(t) (p + V) = (2.21) 
follows. The number of particles per unit physical volume is 

f d 3 P 

n(t)=gJj^f d [a(t)Pf}, (2-22) 



and obeys 

dn(t) 



+ 3H(t)n(t)=0, (2.23) 



dt 

namely, the number of particles per unit comoving volume n(t) a 3 (t) is conserved. 

These are generic results for the kinetic energy momentum tensor and the particle density for any distribution 
function that obeys the collisionlcss Boltzmann equation (|2 . 1 0|) . 
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The entropy density for an arbitrary distribution function for particles that decoupled in or out of LTE is 

d 3 Pf 



s d {t) = 



(2tt) 3 



/ d ln/ d ±(lT/d)ln(lT/d) 



where the upper and lower signs refer to Fermions and Bosons respectively. Since dfd/dt = it follows that 

ds d {t) 



dt 



3H(t) s d (t)=0, 



therefore the entropy per comoving volume Sd{t) a 3 (t) is constant. In particular the ratio 



Y = 



s d (t) 



(2.24) 



(2.25) 



(2.26) 



is a constant for any distribution function that obeys the collisionless Liouville equation [33 |. 

In the case of LTE, using the distribution eq. (|2.14p in the entropy density eq. (|2.24[) yields the result 

Pd + Pd 



s d {t) = 



I'd 



T d a 3 it) T d 



n(t) 



(2.27) 



for either statistics, where pd] P d are evaluated at the decoupling time td- The entropy of the gas of decoupled particles 
does not affect the relationship between the photon temperature and the temperature of ultrarelativistic particles 
that decouple later which can be seen as follows. 

Consider several species of particles, one of which decouples at an earlier time in or out of equilibrium with 
the distribution function f d and entropy given by eq. (|2.24[) while the others remain in LTE with entropy density 
(2 7r 2 /45) g(T) T 3 , until some of them decouple later while ultrarelativistic. Here T is the temperature at time t and 
g(T) is the effective number of ultrarelativistic degrees of freedom. Entropy conservation leads to the relation, 



45 



g(T) T 3 



&d 



a 3 (t) — constant 



(2.28) 



however, because s d (t) a 3 (t) — constant, the usual relation g(T) T 3 a 3 (t) — constant, relating the temperature T of 
a gas of ultrarelativistic decoupled particles to the photon temperature follows. 

For light particles that decouple in LTE at temperature T d ^> m we can approximate 

^/m 2 + p 2 c - fi d p c - [i d _ P f - fi d (t) 



where 



T d 



Td(t) 



Td 



Td(t) 



T d 



Hd{t) 



a(t) ' ^ w a(t) 



(2.29) 



(2.30) 



are the decoupling temperature and chemical potential red-shifted by the expansion, therefore for particles that 
decouple in LTE with T d 3> m we can approximate 



fd(Pf,t) 



1 



1 



1 



e T iW ± 1 e T d<" T i ± 1 



e~^T- ± 1 



(2.31) 



This distribution function is the same as that of a massless particle in LTE which is also a solution of the Liouville 
equation, or collisionless Boltzmann equation. 

Since the distribution function is dimcnsionless, without loss of generality we can always write for a particle that 
decoupled in or out of LTE 



fdiPc 



. I Pc m 

id [ ^r; j^' a ' 



(2.32) 



where are dimensionless constants determined by the microphysics, for example dimensionless couplings or ratios 
between T d and particle physics scales or in equilibrium [i d jT d etc. To simplify notation in what follows we will not 
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include explicitly the set of dimensionless constants m/Td, on, etc, in the argument of fd, but these are implicit in 
generic distribution functions. If the particle decouples when it is ultrarelativistic, m/Td ~ * 0. 



It is convenient to introduce the dimensionless ratios 

Pc _ Pf 



y 



T d T d {t) 



a(t) 



and 



m m 

Xd = Y d ' x ^' = Tdjtj = ' Xd ' 



(2.33) 



(2.34) 



For example, for a particle that decouples in equilibrium while being non-relativistic, fd is the Maxwell-Boltzmann 
distribution function 1321 



fd{Pc) 



2 2 7T2 



45 



Pj 22 TT 2 

3d Yoo e 2 m T " = — 9d Yoc e 
45 



(2.35) 



where gd is the effective number of ultrarelativistic degrees of freedom at decoupling, Y — n/s and is the solution 
of the Boltzmann equation, whose dependence on Xd = m/Td and the annihilation cross section is given in chapter 

5.2 in ref. 

Changing the integration variable in eas. (12. 18|) - i|2.22[) to Pf = y Td{t) we find 



If 00 y 2 
P = 9 ™>T d (t) I p [x] , I p\- x ^=2^J V 2 \j 1 + ^2 fd(y) dy 

V = g^-I v [x] , l v [x] = r dy 4= 
3to 2tt z J q L 



y fd(y) 



3 dl p [x] 



2 7T Z 



leading to the equation of state: 



V 2 fd(y) dy = g Tj{t) I p [x = oo] 
[a;] 1 d In J p [x] 



r i ^ 
»N = - = o 2 r r l 



3 dlnx 

In the ultrarelativistic and non-relativistic limits, x — > and x — > oo, respectively, we find 



1 r r dy . M ^o r 00 y J , 

x Jo ^ fdiv) ' /PN = X J -2^ fd 



(y) 



°° V 2 dy , / x , r , s-oo Z" 00 y 4 dy 



2tt 2 



/ d (y) , Jp[x] = 



In the ultrarelativistic limit the energy density and pressure become, 

r0 ° y 3 ^ 



P = 9 T d (t) 



o 



2tt 2 



Mi/) . p 



P 



2tt 2 



u> x 



/-(!/) 



describing radiation behaviour. In the non-relativistic limit 

P X =°° m g T|(t) 
and the equation of state becomes 



y 2 fd(y)^ =mn(t) , 7> 



3 m 



y 4 fd(y) 



dy 
2tt 2 



w\x\ 



x — >oo 1 



Td{t) 



2 Jp^y 4 dy f d (y) 

/ °° y 2 dy fd(y) 



o, 



corresponding to cold matter behaviour. In the non-relativistic limit, it is convenient to write 



p = mn 7 (t) g 



Td(t) 
T 7 (t) 



Jo y 2 fd,a(y) dy 
4C(3) 



771 77 7 (£) 



Jo y 2 Aa(y) dy 
2<?d C(3) 



(2.36) 



(2.37) 



(2.38) 



(2.39) 



(2.40) 



(2.41) 



(2.42) 
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ga is the number of ultrarelativistic degrees of freedom at decoupling, and n 7 (i) is the 



where C(3) = 1.2020569. 
photon number. 

The average squared velocity of the particle is given in the non-relativistic limit by 



(V 2 



P 2 \ _ Imfc 3 fdHt)Pf] 

2 , 



1)1 



Therefore, the equation of state in thermal equilibrium is given by 



V = ^(V 2 ) P = a 2 p 



a{z) 

and we used that 



Ut) 



Jo y A fd(y) dy 

3 /o°° y 2 fd(y) d v 



0.05124 



1 + z 



9d 



J™y 4 f d (y)dy 
f™y 2 f d (y)dy 



where a is the one dimensional velocity dispersion given at redshift z by 



/keV\ 
V m / 



Jo y 4 fd(y) dy 
/ °° y 2 fd(y) dy 



(t) 



(2.43) 



(2.44) 



(2.45) 



T d {t) = T d (1 + z) = — T 7 (1 + z) 



(2.46) 



T 7 = 0.2348 x 10 _3 eV is the photon temperature today [491 ] . 

The results above, eas. (|2.36|l - f)2.45p are general for any distribution of decoupled particles whether or not the 
particles decoupled in equilibrium. 

Using the relation (|2.42p for a given species (a) of particles with g a degrees of freedom, their relic abundance today 
is given by 



n a h 2 = 



m a g a J °°y 2 fd.aiy) dy 



25.67 eV 



2 9d,a C(3) 



where we used that today h 2 n 7 /p c = l/25.67eV |Hj]. 

If this decoupled species contributes a fraction v a to dark matter, with Q a 
[4^ | for non-baryonic dark matter, then: 



v a &dm and using that il 



DM 



m a 9a / °° y 2 fdAv) dy 



2.695 eV 2 5d , a C(3) 



Since < v a < 1 we find the constraint 



< 2.695 eV 



2 9d,a C(3) 



9a J °°y 2 fdAv) dy ' 



(2.47) 



0.105 



(2.48) 



(2.49) 



where in general fd depends on the mass of the particle as in eq. (|2.32|l . For a particle that decouples while non- 
relativistic with the distribution function eq. (|2.35p this is recognized as the generalization of the Lee- Weinberg lower 
bound [3l|, whereas if the particle decouples in or out of LTE when it is ultrarelativistic, in which case fd,a(y) does 
not dep end on the mass, eq. (|2.49[) provides and upper bound which is a generalization of the Cowsik-McClelland 
[H Hof bound. 

The constraint eq. (|2.49p suggests two ways to allow for more massive particles: by increasing gd, namely the 
particle decouples earlier, at higher temperatures when the effective number of ultrarelativistic species is larger, and/or 
decoupling out of LTE with a distribution function that favors smaller momenta, thereby making the denominator 
in eq. (!2.49|) smaller, the smaller number of particles allows a larger mass to saturate the DM abundance. 

For the particle to be a suitable dark matter candidate, the free streaming length must be much smaller than the 
Hubble radius. Although we postpone to a companion article [5l[ a more detailed study of the free streaming lengths 
in terms of the generalized distribution functions, here we adopt the simple requirement that the velocity dispersion 
be small, namely the particle must be non-relativistic 



{v 2 ) = A)«i. 



(2.50) 
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From ea. (|2.43|) this constraint yields 



m ^ U^v i fd{y)dy 
t&) V/0 00 v 2 fdv)dy ' ( j 

where T d (t) is given by eg. ([235]) . From eqs. ([2TT5]) . piM]) . f27jj]) and (|2~5Tj) we obtain the following condition for the 
particle to be non-relativistic at redshift z 



to > 2.958 — — x 10 4 4/ 700 eV . (2.52) 

ff | V Jo 2/ Jd(y) e?2/ 

Taking the relevant value of the redshift for large scale structure to be the redshift at which reionization occurs z s ~ 10 
[52| . we find the following generalized constraint on the mass of the particle of species (a) which is a dark matter 
component 



2-958 4 Jo y*fd(y)dy 2g d ((3) 

— — x 10 «/ -755 9 , 7 r , eV < to < 2.695 x eV . (2.53) 

g* V Jo V fd{y) dy gj Q y 2 fd(y)dy 

The left side of the inequality corresponds to the requirement that the particle be non-relativistic at reionization 
(taking z s ~ 10), namely a small velocity dispersion (V 2 /c 2 ) <C 1, corresponding to a free streaming length A/ s ~ 



(V 2 /c 2 ) dn much smaller than the Hubble radius (djj), while the right hand side is the constraint from the 
requirement that the decoupled particle is a dark matter component, namely eq. (|2.49[) is fulfilled. 



III. LIGHT THERMAL RELICS AS DARK MATTER COMPONENTS. 



In this section we consider particles that decouple in LTE. 



A. Fermi-Dirac and non-condensed Bose-Einstein gases of light particles as DM components. 

The functions I p (x), I-p(x) in the density and pressure denoted by I±(x), J±(x) respectively for Fermions (+) and 
Bosons (— ) and the equation of state w[x) eq. (|2.37p for each case are depicted in figs. [TJH for vanishing chemical 
potential in both cases. We have also numerically studied these functions for values of the chemical potential in the 
range < \fid\/Td < 0.5 but the difference with the case of vanishing chemical potentials is less than ~ 5% even for the 
largest value studied \nd\/Td = 1.0 which is about the maximum consistent with constraints on lepton asymmetries 
allowed by BBN and CMB [H. 




FIG. 1: Fermions without chemical potential. Left panel: I+(x) and J+(x) vs x. Right panel: 3 w[x] vs. x. 1+ = I p , J+ = I-p- 
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FIG. 2: Bosons without chemical potential. Left panel: I-(x) and J-(x) vs a;. Right panel: 3 w[x] vs. x. 7_ = I p , J- = I-p. 



These figures make clear that the onset of the non-relativistic behavior occurs for X'n/f" ^ 5 111 both cases. It is useful 
to compare this result, with the generalized constraint eq. (|2.51[) for the case of thermal relics. Replacing the LTE 
distribution functions (Fermi-Dirac or Bosc-Einstcin, without chemical potentials) in cq. (|2.51[) we obtain 



> 3.597 for Fermions 



> 3.217 for Bosons 



(3.1) 

5 in both cases for the 



The detailed analysis of the corresponding functions yields the more precise estimate x nr > 
transition to the non-relativistic regime. 

Therefore, the decoupled particle of mass m becomes non-relativistic at a time t* when m > 5T^(t*). At the time of 
Big Bang Nucleosynthesis (BBN) when [32[ Tbbn ~ O.lMeV and gsBN ~ 10, the decoupled particle is non-relativistic 
if 



> 



g d 3 MeV , 



(3.2) 



in which case it does not contribute to the effective number of ultrarelativistic degrees of freedom during BBN and 
would not affect the primordial abundances of light elements. If the particle remains ultrarelativistic during BBN the 
total energy density in radiation is 32] 



ptot(t) = gQ g*(t) T*{t) 



eg ( g*(t) 



9*(t) 



0,1 



(3.3) 



where T 7 is the (LTE) temperature of the fluid, c = 1(7/8) for Bosons (Fermions), g*(t) is the effective number of 
ultrarelativistic degrees of freedom at time t from particles that remain in LTE at this time, and gd is the effective 
number of degrees of freedom at decoupling. The second term in eq. (|3.3p is an extra contribution to the effective 
number of ultrarelativistic degrees of freedom. 

At the time of BBN, g*{tBBN) ~ 10 [13] and early decoupling of the light particle, gd ^> g*(tBBN), leads to a 
negligible contribution to the effective number of ultrarelativistic degrees of freedom well within the current bounds 
[541 ]. Therefore, provided that the decoupled particle is stable, either for light particles that remain relativistic during 
(BBN) but that decouple very early on when gd 3> 10 or when the particle's mass m > lMeV, there is no influence 
on the primordial abundance of light elements and BBN does not provide any tight constraints on the particle's mass. 



B. A Bose condensed light particle as a Dark Matter component 

Consider the case of a light bosonic particle, for example an axion-like-particle. Typical interactions involve two 
types of processes, inelastic reactions are number-changing processes and contribute to chemical equilibration, while 
elastic ones distribute energy and momenta of the intervening particles, these do not change the particle number but 
lead to kinetic equilibration. Consider the case in which chemical freeze out occurs before kinetic freeze-out, such is 
the case for a real scalar field with quartic self-interactions. In this theory, number-conserving processes such as 2 ^ 2 
establish kinetic (thermal) equilibrium, but conserve particle number, a cross section for such process is oc A 2 where 
A is the quartic coupling. The lowest order number-changing processes that contribute to chemical equilibrium are 
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4 ±5 2, with cross sections oc A 4 . Hence this is an example of a theory in which chemical freeze out occurs well before 
kinetic freeze out for small coupling. 

Another relevant example is the case of WIMPs studied in ref. [33[ where it was found that T c d ~ 10 GeV, while 
Tkd ~ lOMeV where T c d,Tkd are the chemical and kinetic (thermal) decoupling temperatures respectively. Although 
this study focused on a fermionic particle, it is certainly possible that a similar situation, namely chemical freeze-out 
much earlier than kinetic freeze out, may arise for bosonic DM candidates. 

Under this circumstance, the number of particles is conserved if the particle is stable, but the temperature continues 
to redshift by the cosmological expansion, therefore the gas of Bosonic particles cools at constant comoving particle 
number. This situation must eventually lead to Bose Einstein condensation (BEC) since the thermal distribution 
function can no longer accomodate the particles with non-vanishing momentum within a thermal distribution. Once 
thermal freeze out occurs, the frozen distribution must feature a homogeneous condensate and the number of particles 
for zero momentum becomes macroscopically large. Although some aspects of Bose Einstein condensates were studied 
in ref. (35l . |3(| , we study new aspects such as the impact of the BEC upon the bound for the mass and the velocity 
dispersion of DM candidates. 

The bosonic distribution function for a fixed number of particles includes a chemical potential and is given by 
eq. (l2.14[) where fid < m for the distribution function to be manifestly positive for all p. Separating explicitly the 
contribution from the p = mode the number of particles per comoving volume V c is 



1 



1 



- 1 



1 



n 



d 3 p c 



1 



(2tt) 3 y m 2+pg-^ 

e T d - 1 



(3.4) 



where 



n 



1 1 

Vc 



e T d - 1 



(3.5) 



is the comoving condensate density. In the infinite volume limit the condensate term vanishes unless fid — > m. For 
m/Td < 1 we find 



n = n + 



n c(3) 



' HA 



(3.6) 



where 



ltd. 

e T d 



C(3)V *' 



(3.7) 



The maximum value that fid can achieve is to, therefore, neglecting m/Td we replace Z[e T d] by Z[l] = 1. If the 
comoving particle density 



n > 



ric(3) 



(3.8) 



then, there must be a zero momentum condensate with no ^ and fid = m in the infinite (comoving) volume limit. 
In this limit we find, 



n { 



for T d <T c 
for T d >T c 



(3.9) 



where the critical temperature is given by 



C(3) 



(3.10) 



The solution of the equation (|3.4p that determines the condensate fraction shows that for Td < T c 



fi d = m 



(3.H) 
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In the infinite volume limit the distribution function for particles that decouple while ultrarelativistic m/T c i <C 1, for 
Td < T c becomes 



fd{p c ) = n 5 (3) (p c ) + Ti^ ■ 

e T d - 1 

From eq. ()2.22p the total number of particles for m/Td <1, T c > Td is given by 

n(t)=n (t) + ^-Tl(t), 



(3.12) 



(3.13) 



where 



For T d < T c eq.pJ]) implies that 



n (t) 



no 



^] -1 



hence for < T c the total density is given by 



(3.14) 



(3.15) 



(3.16) 



The enhancement factor (T c /Td) 3 over the thermal result reflects the population of particles in the condensed, zero 
momentum state. The energy density and pressure are given by 



p{t) = gm{n Q (t)+T 3 (t)I^[x(t)}} 
T 5 (t) 

V{t) = g IWW)] 
6 m 



where 



i 

2^2 

1 

2^2 



I ^ J- ^— dy 



x 2 {t) ev - 1 



a/1 



— dy , x(t) 



■ 



(3.17) 
(3.18) 

(3.19) 
(3.20) 



are the contributions from the particles outside the condensate (p ^ 0). 

Two important aspects emerge from these expressions: i) the condensate always contributes as a non-relativistic 
component, ii) the condensate does not contribute to the pressure. 

Replacing eq. (|3.15p into (|3. 1T[) and using 



V 2 dy 



= 2C(3), 



ey - 1 

the energy density and equation of state for Td < T c can be written compactly as 

p(t) = gm T[x{t)} 

where 

3 



'1 



V 



y 



x 2 {t) I ev - 1 



dy 



(3.21) 



(3.22) 



The equation of state 



w\x\ 



3x 2 l[x] 



(3.23) 
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is displayed in fig. [31 from which it is clear that for T c /Td > 1 the non-relativistic limit sets in much earlier than for 
the non-Bose condensed case. This is a consequence of the zero momentum particles in the BEC which contribute as 
pressureless cold matter, even when the light bosonic particle decouples while ultrarelativistic. 

For Td < T c when the particle becomes non-relativistic, namely x — > oo, the energy density becomes 

p(t) = g m n{t) (3.24) 

where n(t) is the total number of particles per physical volume, including the condensate and non-condensate com- 
ponents, from eq. (|3.16p it follows that 

p(t)=gm^(^j 3 T*(t), (3.25) 
from which for T c > Td it follows analogously to eq. (|2.47|) that 



n BE h 2 = -(—) . (3.26) 

25.67 eV g d \T d J 



3 



The dark matter fraction that these particles can contribute is given by 

m g_ ( TA 
2.695 eV g d \T d ) 



^ = ^- T ^m\ (3.27) 



resulting in the upper bound 



9d (T d ^ 3 



m < 2.695 — -f ) eV . (3.28) 
9 V-'c 



In the Bose condensed case Td/T c < 1 the bosonic particle is light unless it decouples very early on at high temperature 
with a large g d - The presence of a BEC tightens the constraint on the mass of the light bosonic particle via the extra 
factor (T d /T c f in fsjZS]) . 

Of importance for clustering, is the velocity dispersion when the particle becomes non-relativistic, it is given by 



i 2 / C(3) 



T d (t) 



(3-29) 



where C(5) = 1.0369278.... 

The presence of the BEC, accounted for by the factor (^Td/T^j < 1 in cq. (|3.29[) , diminishes the velocity dispersion. 

This is a consequence of the fact that the particles in the condensate all have vanishing momentum, and only the 
non-condensate particles contribute to the velocity dispersion but the fraction of particles outside of the condensate 
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is precisely the factor [T d /T c j . Therefore the presence of a BEC leads to a decrease in the velocity dispersion and 
consequently even for light particles to a decrease in the free streaming length. 

These results imply that although the Bosonic particle is bound to be very light by the bound (|3.28[) (unless they 
decoupled very early ), if T c 3> T<j it is not a HDM component but can effectively act as either a WDM or CDM 
because of a small velocity dispersion. Whether T c 3> T d or not has to be studied within the microscopic particle 
physics model that describes this DM component. 



IV. COARSE GRAINED PHASE SPACE DENSITIES AND NEW DM BOUNDS 

In their seminal article Tremaine and Gunn [3(| argued that the coarse grained phase space density is always 
smaller than or equal to the maximum of the (fine grained) microscopic phase space density, which is the distribution 
function. Such argument relies on the theorem [43| that states that collisionless phase mixing or violenty relaxation by 
gravitational dynamics can only diminish the coarse grained phase space density. A similar argument was presented 
by Dalcanton and Hogan [Til . [3q , and confirmed by recent numerical studies [Ij] . 

As noticed in ref . [361 ] , the case of the Bose-Einstein distribution, requires a careful treatment because for massless 
particles the Bose-Einstein distribution diverges at small momentum. This divergence is present if there is a BEC 
even when the mass of the Bosonic particle is included. This is so since p d = m is required to form a BEC and 
the distribution functions diverge at zero momentum, even the part of the distribution function that describes the 
particles outside the condensate diverges at P — 0. Madsen recognized this caveat in the Bosonic case and in ref. [36| 
introduced an alternative statistical interpretation of the phase space density, similar to that introduced in [Til . |38| but 
with the upper limit in the momentum integrals replaced by a (physical) momentum cutoff as suggested by the phase 
mixing theorem 

El- 

However, it is straightforward to show that the resulting coarse grained phase space density is 
not a Liouville invariant. Instead, we define the coarse grained (dimcnsionlcss) primordial phase space density 



V 



(4.1) 



which is Liouville invariant and where (-P 2 ) is defined in cq. (|2.43p . Since the distribution function is frozen and is a 
solution of the collisionless Boltzmann (Liouville) equation (|2.1ip it is clear that I? is a constant, namely a Liouville 
invariant in absence of self-gravity, . Including explicitly a possible BEC, T> is given by 



V = 



2tt 2 



27T Tip 



+ XT y 2 fd(y)dy 



(4.2) 



where fd(y) is the distribution function for the non-condensed particles in the Bosonic case and no is the comoving 
density of the Bose-Einstein condensate 



2?r 2 n _ J 2C(3) 




T 3 
± d 



i 



for the BEC with T d < T c 



(4.3) 



for the Fermionic or non — Bose condensed case 



When the particle becomes non-relativistic p(t) = m n(t) and (V^ 2 ) = (— therefore, 

p Qdh 



V 



(4.4) 



where Q dh is the phase-space density introduced by Dalcanton and Hogan [ll|, [36 

P 



Q 



DH 



(y 2 Y 



and the one-dimensional velocity dispersion a is defined by eq. (|2.44|) . 



(4.5) 
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In the non-relativistic regime V is related to the coarse grained phase space density Qtg introduced by Tremaine 
and Gunn (30| 

3 

Qtg = t J . = ( £Y * ■ ( 4 - 6 ) 



to 4 (2 7r(7 2 )! 



2tt 



The observationally accessible quantity is the phase space density p/a 3 , therefore, using p = m n for a decoupled 
particle that is non-relativistic today and ea. (|2.44|) . we define the primordial phase space density 



2 

PDM - 3* m 4 V = 6.611 x 10 s V 



a 3 
"dm 



keV 



4 Mo/kpc 
(km/: 



s 



(4-7) 



where we used that keV 4 (km/s) 3 = 1.2723 10 8 

During collisionless gravitational dynamics, phase mixing increases the density and velocity dispersions in such a 
way that the coarse grained phase space density either remains constant or diminishes, namely 



A; < 6.611 x 10 8 V 
a A 



keV 



M ° /kpC 3 3 . (4.8) 
(km/s) 



where T> is given by eq. (|4.2p for an arbitrary distribution function. For a particle that decouples when it is ultrarela- 
tivistic T> does not depend on the mass, hence eq. (|4.8p yields a lower bound on the mass of the particle directly from 
the observed phase space density and the knowledge of the distribution function. 

For comparison it is convenient to gather the values T> eq. (|4.2p for the usual LTE cases that follow from eqs. (|2.31l) 
and [g!51l 

1.963 x 1CT 3 Fermions, fi d = 

3.657 x 1CT 3 Bosons without BEC 

V = gx { ,/ T ^ (4.9) 

' 3.657 x lO- 3 f^-J Bosons with BEC, T c > Td 

8.442 x 10~ 2 gd Yoo non — relativistic Maxwell — Boltzmann. 

where gd is the number of ultrarelativistic degrees of freedom at decoupling. 

We note that for T c 3> Td the presence of a BEC increases dramatically the primordial phase space density. This is 
a consequence of the enhancement of the particle density over the thermal case due to the presence of the condensate, 
and the decrease in the velocity dispersion because the particles in the condensate all have zero momentum. 



A. New bounds from phase space density and dShps-data 

We derive here new bounds from the latest compilation presented in ref. [0j directly on p/a 3 for the dataset 
comprising ten satellite galaxies in the Milky- Way dSphs. It proves convenient to write eq. (|4.8[) as 



m 4 > 



62.36 eV 



V 



1Q _4 P_ ( km /s) 

T 3 m /i,„„3 



a 3 M Q /kpc J 



(4.10) 



the data in ref. [lj] yields the range 



!M> 1(T 4 A; ( km//S ) „ 20 . 

t-3 n/r n A — 1 



ct 3 M Q /kpc 3 



(4.11) 



and we choose a fiducial value for this quantity in the middle of the range of the data [l4[ ~ 5 — 10, leading to the 
new bound 



m > — r ev . 



(4.12) 
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For thermal relics that decoupled while ultrarelativistic with vanishing chemical potentials and no BEC, we find from 
eqs. flU} and (|332]1 . 



> 



1 I 0.475 keV Fermions 
A 1 0.407 keV Bosons without BEC 



(4.13) 



and for Bosons with BEC (T d < T c ) we find 



to > -4 0.407 keV 

S 4 



Bosons with BEC 



(4.14) 



For particles that decouple out of LTE with arbitrary distribution functions the form of the new bound is given by 
eq. (|4.12p with V given by eq. (|4.2|) . The detailed form of V is completely determined by the distribution function at 
decoupling, which must be obtained from a microscopic calculation of the kinetics of decoupling. Once the distribution 
function is obtained, the new bound cq. (|4.12[) yields the lower bound of the mass consistent with the observational 
data. 

Combining the upper bound (|2.49[) with the lower bound eq. (|4.10p we establish the mass range for the DM candidate 



62.36 eV 



10" 



(km/s) 
M /kpc 3 



< to < 2.695 eV 



2 g d C(3) 



9 I^v 2 fd(y) d v ' 



(4.15) 



where T> is given by eq. (|4.2|l and the compilation of data in [14( constrains the bracket [• • • ] 4 ~ 1 — 2. 

For thermal relics that decoupled in LTE while ultrarelativistic, and taking the bracket in the middle of the range 
we obtain from eqs. (|4.9[) and (14.15[) 



444 eV 






9d 


1 


< 


TO < 


S 4 






9 


380 eV 






9d 


1 


< 


TO < 


9 4 






9 



380 eV 



9 3 



Td 

T r 



< to < — 2.695 
9 



III 



I 3 



eV BEC 



(4.16) 



Therefore, if the thermal relic decouples in equilibrium this mass range indicates that it must decouple when gd g 4 > 
110 — 150, namely at or above the electroweak scale [32j. In the BEC case, for Td <C T c the fulfillment of the bound 
requires very large g d 5~ 4 , namely thermal decoupling at a scale much larger than the electroweak scale. 

An alternative is that the particle is very weakly coupled to the plasma and decouples away from equilibrium with 
a distribution function that yields a smaller abundance increasing the right hand side of eqn. (|4.16l ). 



B. Generalized Tremaine-Gunn bound. 



The Tremaine-Gunn bound [30( establishes a relation between the properties of dark matter in galaxies through 
their phase space densities. It assumes that dark matter could be reliably described by an isothermal sphere solution 
of the Lane-Emden equation with the equation of state (|2.44p [HH [5^ | . In thermal equilibrium the quantity [56| 



n 



G m 2 N 2 G p L 2 



L T 



3 a 2 



(4.17) 



is bound to be i] < 1.6 to prevent the gravitational collapse of the gas. Here V = I? stands for the volume occupied 



by the gas, N for the number of particles and T 



to a for the gas temperature. The length L is similar to 



the King radius [55j. However, the King radius follows from the singular isothermal sphere solution while L is the 
characteristic size of a stable isothermal sphere solution (56[. 

Combining eq. (|4.17[) with eq. (|4.8|) results in a generalized Tremaine-Gunn bound 



4 \ 
TO > 



[85.22 eV] 4 10 km/s 



2 G L 2 a V 



V 



kpc 



(4.18) 
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therefore the generalized Tremaine-Gunn bound on the mass becomes 



85.22 eV i 
m > i r\ A 



2?4 



10km/s 


4 


kpc 


a 




L 



(4.19) 



The compilation of recent photometric and kinematic data from ten Milky Way dSphs satellites [T3| yield values for 
the one dimensional velocity dispersion (<r) and the radius (L) in the ranges 



0.5 kpc < L < 1.8 kpc , 6.6 km/s < a < ll.lkm/s . 



(4.20) 



For particles that decouple in LTE when they are ultrarelativistic (ultrarelativistic thermal relics) with vanishing 
chemical potential and no BEC we find from eas. (|4.9p and (|4.19|) . 



m > I — 

9 



10 km/s 


4 


kpc 




a 




L 





0.405 keV Fermions 
0.347 keV Bosons . 



For the case of ultrarelativistic bosonic thermal relics with a BEC and T d < T c we find the bound 
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m > 



0.347 keV 



10 km/s 


4 


kpc 


2 


~T d 


a 




L 




T c 



(4.21) 



(4.22) 



Therefore, the BEC case allows for smaller masses to saturate the Tremaine-Gunn bound for T c 3> T d , a consequence 
of the enhanced primordial phase space density in the presence of the BEC. 



C. DM mass values from velocity dispersion 



We can use the independent data provided in ref. 14] on the mean density and velocity dispersion to explore 
bounds solely from the velocity dispersion. Since the phase space density only diminishes or remains constant during 
the collisionless gravitational dynamics of clustering, from which it follows that 



PDM Ps 

a dm ~ a l 



(4.23) 



where pdm and o~dm are, respectively, the matter density and velocity dispersion of the homogeneous dark matter 
prior to gravitational collapse. p s and o~ s are, respectively, the satellite's mean volume mass density and velocity 
dispersion. Assuming that DM has a single component, its density today is [49j | 



Pdm 



DM 



h 2 1.054 x 10 4 eV/cm 3 = 1.107 x 10 3 eV/cm 3 , 



(4.24) 



o~dm is given by eq. (|2.45p . Ref. [14| quotes the following values for the favored satellite's cored dark matter density 
and velocity dispersion 



GeV 



10 



km 



(4.25) 



Eqs. (|4~2Uj) . (|4~M|) and (|4~2"5)) lead to 



o-dm < 0.06 



km 



(4.26) 



Combining eq. (|2.45|) for 



and eq. p~2"6]) yields 

0.847 



keV 



> 



Jo°° y 4 My) dy 
Jo°° y 2 fd(y) d y 



(4.27) 



For thermal fermions or bosons without chemical potential (no BEC) and 10 < g& < 100 we find m ~ 0.6 — 1.5 keV 
in agreement with the bounds found above and the conclusions of ref. [57j | . A suppression factor (Td/T c ) 3 appears in 
the BEC case for the same range of g d - 
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We emphasize that the bound eq. (|4.27p is independent from the bound eq. (l4.13p obtained from the phase space 
density above, and relies on the fact that the observational data [hij yields separate information on p s and a s . 

It proves illuminating to analyze the velocity dispersion udm from expression (|2.45| at z — for thermal relics. 
We find 

! 0.187 Fermions pd = 
0.167 Bosons no BEC 
0.167 x Bosons with BEC, T c > T d 

0.09 x Jx2 non — relativistic 



We see that for T c ^S> Td light Bosonic particles that decoupled while ultrarelativistic but undergo BEC can effectively 
act as CDM with very small velocity dispersion. 

In ref. [33| it is found that kinetic decoupling for a WIMP of mass m ~ 100 GeV occurs at Td ~ 10 MeV, leading to 
the estimate ^fxd = \JmjTd ~ 100. Thus, for CDM from weakly interacting massive particles the velocity dispersion 
cq. (l4~28l) is: 

a /100GcV\ ^fxd -i / km\ 
<r wimp ~ 10- (_^-) ^ 9 /(9 7 ). (4.29) 

Thus, (Jwimp is eight orders of magnitude smaller than the typical velocity dispersion in dSphs [li[ for wimps of 
m ~ 100 GeV that decoupled in LTE at T d ~ 10 MeV [Hj]. 

It is noteworthy to compare the phase space densities of the homogeneous dark matter distribution for the thermal 
relics that decoupled ultrarelativistically and non-relativistically with that observed in the satellites dShps. If the 
distribution of dark matter is cored [lj] 1 

^) -sxio 6 ^/^ . (430) 



I km/s 



If the distribution of dark matter is cusped, ref. [14J gives the value for the density p s ~ 2 TeV/cur yielding 

~2xlO»^L . (4. 31) 



^ s / cusped 



Assuming that a thermal relic that decoupled when ultrarelativistic is the only DM component with the density given 
by the value today p DM ~ 1-107 x 10 3 eV/cm 3 HU, we find from eqs.^MB at z = 0, 

v/ 3 / \ 3 ( 0.177 Fermions 
BEE. ^ i 6 ^ V/Cm 3 J!L j gd ) 0.247 Bosons without BEC (4.32) 
'km/s] \ kcV / I 0.247 (T c /T d ) I Bosons with BEC . 



Thus, for g d > 10 we see that for m ~ keV the phase space density for thermal relics that decoupled being ultrarel- 
ativistic is of the same order as the phase space density in dShps with cores, eq. (|4.30p . Thermal relics with mass in 
the ~ keV range obviously favor cores over cusps because the primordial phase space is Pdm/o~i) M > p s /c 3 for cores 
while pdm I >0 ~dm ^ Ps/^s f° r a cuspy distribution, and according to the theorem in [30l. [43|. the phase space density 
can only diminish during gravitational clustering. 

An enhancement factor (T c /Td)^ appears in eq. (|4.32|) for the case of a BEC. Notice, that for T c /T d > 10 and 
m ~ keV, a BEC yields a phase space density consistent with cusps as a result of the small velocity dispersion and 
the CDM behavior. 

Recent TV-body simulations [3] indicate that the phase space density decreases by a factor 10 — 10 2 due to grav- 
itational relaxation during structure formation between < z < 10, with smaller relaxation in WDM than in CDM 
[ill [lil]. Therefore, from these numerical results it follows that 



EL „ 10- 2 . (4.33) 



DM 



1 (eV/c 2 )/cm 3 = 0.026 M /kpc 3 
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Combining this result with the observational results eqs. (|4.30|) - (|4.31[1 and the primordial phase space density eq. (|4.32|) 
for a thermal relic that decoupled while ultrarelativistic, we find 

mcored ~ ~~r keV , m cusp — keV . (4.34) 

These values and the upper bounds for m in eqs. (|4.16p yield the following bounds for thermal relics 

9d <7 3 > 500 for cores , gd g~ z > 2000 for cusps . (4.35) 

Therefore, thermal relics, DM candidates that decouple when relativistic, must decouple at a temperature well above 
the electroweak scale. Eqs. (|4.34|) and (|4.35|) imply for the mass value: 

2 8 
m CO red ~ — keV , m cusp ~ — keV . (4.36) 
9 4 3 4 

Although m cusp is not too much larger than m core d it is noteworthy that the thermal relic DM candidate that leads 
to cusped profiles must decouple when ga > 2000 namely very early at a temperature scale corresponding to a grand 
unified theory with a large symmetry group. 

For the case of CDM from wimps which decoupled non-relativistic, we find from eqs. (l4.24p and (|4.28[) 

Pvimp in24 eV/cm 3 ( m \ ( 100 V 

TZ^y [^Gev) feJ (4 ' 37) 



km/s 



The phase space density always decreases by dynamical relaxation, a result recently confirmed numerically by iV-body 
simulations 44j. For initial values of the phase space density which are much lower than the primordial ones, these 
yield a typical decrease by a factor 10 2 — 10 3 (44J. If these results should persist in N-body simulations with larger 
values of the initial phase space density, they would imply a tension between the phase space density of WIMPs 
eq. (|4.37p being eighteen to fifteen orders of magnitude larger than that in dShps either cored eq. (|4.30|) or cusped 
ea. (l43il r 14|. 

Combining eqs. (|430|) . P~5Tj) , (1435)) and ((4~3Tj) yield for wimps as DM, 

\ZmTd ~ ^ keV for cores , \f mTd ~ keV for cusps . (4.38) 

9l gJ 

^,From the combined analysis of the primordial phase space densities, the observational data on dSphs [l4| and the 
iV-body simulations in ref. [4J|j we conclude the following: 

• (i): Thermal relics with m ~ few keV that decouple when ultrarelativistic lead to a primordial phase space 
density of the same order of magnitude as in cored dShps and disfavor cusped satellites for which the data [lij 
yields a much larger phase space density. 

• (ii) : Light bosonic particles decoupled while ultrarelativistic and which form a BEC lead to phase space densities 
consistent with cores and if T c /T,i > 10, also consistent with cusps. However, for thermal relics to satisfy the 
bound eq. (|4.16p they must decouple when g^ g~* (Td/T c )« > 130, namely above the electroweak scale. Recall 
that typically g takes a value between one and four. 



V. NON-EQUILIBRIUM EFFECTS: 



The main results of our analysis are the new bounds from DM abundance and phase space density of dShps 
summarized in eq. (|4.15p . When the dark matter candidate decouples out of LTE these bounds establish a direct 
connection with the microphysics via the frozen distribution functions. These functions must be obtained from a 
detailed calculation of the microscopic processes that describe the production and pathway towards equilibration of 
the corresponding dark matter candidate. If kinetic (and chemical) freeze out occur out of LTE the distribution 
functions will keep memory of the initial state and the detail of the processes that established it. 

Non-equilibrium effects have been mainly considered for massive particles that decoupled when non-relativistic [58[ 
or as distortions in the neutrino distribution functions during BBN |59l.[60|. Instead, we focus here on DM constraints 
from decoupling out of (LTE) at temperatures larger than the BBN scale and when particles are ultrarelativistic. 
Decoupling out of LTE in this case has been much less studied. In this section we explore a cosmologically relevant 
mechanism of production and equilibration which describes a wide variety of situations out of LTE. 
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A. Particle production followed by an UV cascade: 



Early studies of particle production via parametric amplification and oscillations of inflaton-like scalar fields revealed 
that particles are produced via this mechanism primarily in a low momentum band of wave vectors (45j leading to a 
non-thermal spectrum (figs. 2-3 in ref. [45| illustrate these effects). 

Subsequent studies [46| showed that the early phase of parametric amplification and particle production is followed 
by a long stage of mode mixing and scattering that redistributes the particles: the larger momentum modes are 
populated by a cascade whose front moves towards the ultraviolet akin to a direct cascade in turbulence, leaving in 
its wake a state of nearly LTE but with a lower temperature than that of equilibrium [46| . 

The dynamics during the cascade process diminishes the amplitude of the distribution function at lower momenta 
and populates the higher momentum modes. The distribution function develops a front that moves towards the 
ultraviolet. Behind the front the distribution function is nearly that of LTE with a different temperature and amplitude 
and slowly evolves towards thermal equilibrium [46| . If these particles are very weakly coupled to the plasma it is 
possible that the advance of the cascade and the front of the distribution towards larger momenta is interrupted when 
the rate of scattering or mode mixing becomes smaller than the expansion rate. In this case, the distribution function 
is frozen well before reaching complete LTE resulting in a population of modes primarily at lower momenta up to the 
scale of the front. This study suggests the following frozen distribution function 



fd{y) = fa f 



y), 



(5.1) 



where feq(-^) is the equilibrium distribution function for an ultrarelativistic particle at an effective temperature 
£ Td- Namely, £ = 1 at thermal equilibrium and £ < 1 before thermodynamical equilibrium is attained. 

This form describes fairly accurately the cascade with a front that moves towards the ultraviolet, which is interrupted 
at a fixed value of the momentum, identified here to be p° c = y T^; Td is the temperature of the environmental degrees 
of freedom that are in LTE at the time of decoupling. 

The amplitude fo and effective temperature £ Td < Td reflect an incomplete thcrmalization behind the front of 
the cascade and determine the average number of particles in its wake (46} . This interpretation is borne out by the 
detailed numerical studies in ref. [46| . For Fermi-Dirac ultrarelativistic particles (with vanishing chemical potential) 
< fo < 2 whereas for Bose-Einstein ultrarelativistic particles < fo < oo. Neglecting the possibility of a BEC, for 
a fermionic or bosonic equilibrium distribution function / eg , we find 
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and the primordial phase space density becomes 
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(5.3) 



where V eq is the phase space density eq. (|4.2j) for the equilibrium distribution /, 



I eq- 



For a fermionic species without chemical potential [f eq {y) = l/(e y + 1)], the bound eq. (|4.15|) becomes 
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and the one dimensional velocity dispersion eq. (|2.45p becomes today: 
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(5.5) 



The functions F(s), H(s) and v(s) for the case f eq (y) = l/(e v + 1) ar e displayed in fig. 2J For the Bose-Einstein case 
without a BEC the behaviors of F(s), H(s) and v(s) are qualitatively similar. 

It is clear that the bound eq. (|5.4p for the range of m can easily be satisfied for moderate values gd ~ 10 — 50 
corresponding to decoupling temperatures lMeV < T4 < 1 GeV and fo £ 3 F(s) < 0.08. 

Remarkably, the non-equilibrium distribution eq. (|5.1D turns out to be a generalization of several non-equilibrium 
distribution functions of cosmological relevance proposed in the literature: 




(a): sterile neutrinos produced non-resonantly via the Dodelson-Widrow mechanism [18l 
bution function is obtained from (|5.ip by taking £ = 1; s — * oo, /o ~ 0.043 keV/m (ref. 
find for the mass range, phase space density and velocity dispersion respectively, 
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(5-6) 

The major uncertainty is the evaluation of g^. In the Dodelson-Widrow [18| scenario the sterile neutrino 
production peaks at T ~ 130 MeV, this temperature is very near the region where the QCD phase transition 
occurs at which the effective number of ultrarelativistic degrees of freedom changes dramatically. // decoupling 
occurs at a temperature higher than the QCD critical temperature, then g& ~ 30 and the mass bound ea. (|5.6[) 
may be fulfilled, but for a lower decoupling temperature when gd < 25 — 30 the mass bound may not be fulfilled. 
If the mass bound is fulfilled, Pdm/^dm ^ s compatible with cored dSphs [3] [see eq. (|4.30|) ] but not with the 
cusped distributions, [see eq. (|4.31|) ]. Combining the bound eq. (|5.6[) . the observed phase space density eq. (|4.30|) 
[lH and the ./V-body results of ref. [3] which yield phase space relaxation by a factor ~ 10 2 we find that 



to ~ — - keV 

9* 



(5.7) 



(b) : sterile neutrinos produced by a net-lepton number driven resonant conversion studied by Shi and Fuller 
[19( for which the distribution function is obtained from eq. (|5.1[) for £ = 1, s ~ 0.7, /o ~ 1, (see fig. 1 in the 
first reference in [l9(). We find 
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Again, a source of uncertainty is the evaluation of g^ because in the resonant-mediated sterile neutrino pro- 
duction, the maximum production rate is near the QCD temperature (l9j . However, it is clear that in this case 
the mass bound is less sensitive to the uncertainty in ga (a small value gd ~ 10 fulfills the bound), although 
the MSW resonance occurs also near the QCD critical temperature [l9j]. The velocity dispersion is small be- 
cause the distribution is skewed towards small momenta. Again, Pdm / <J DM is consistent with cored dSphs [see 
eq. (l4.30p ] but not with cusped distributions [see eq. (|4.31[) ]. A similar analysis as in the previous case combining 
the observational data, the results of ref. [44( and the bound eq. ()5.8() yields 



7)1 



^keV. 

9 i 



(5.9) 



• (c): Our distribution function eq. ()5.ip for £ 
in ref. [H| to model WDM [eq.(8) in ref. [2f 
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The parameter j3 cannot be too small, although a small (3 increases the mass, it decreases the phase space 
density. 

Although recent studies (6lj suggest that the description of the production mechanism of sterile neutrinos must 
be reassessed with likely implications on their distribution functions after decoupling, the above estimates provide a 
guidance to the range of mass, primordial phase space density and velocity dispersions for sterile neutrinos as possible 
WDM candidates. 



VI. CONCLUSIONS 



We have obtained new constraints on light DM candidates that decoupled while ultrarelativistic in or out of LTE in 
terms of their distribution functions. The only assumption is that these distribution functions are homogeneous and 
isotropic. A Liouville invariant coarse grained primordial phase space density is introduced that allows to combine 
phase space density arguments with a recent compilation of photometric and kinematic data on dSphs galaxies to 
yield new constraints on the mass, velocity dispersion and phase space density of DM candidates. The new constraint 
on the mass range is 
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where the primordial phase space density is given by 
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fd{Pc/Td) is the distribution function at decoupling, g the number of internal degrees of freedom of the particle, and 
p/<r 3 is the phase space density obtained from observations. The upper bound arises from requesting that the DM 
candidate has a density < pdm today, and the lower bound arises from requesting that the phase space density in 
halos p/a 3 be smaller than or equal to the primordial phase space density of the collisionless non-relativistic (today) 
DM component 



Pdm/<J 3 D m 
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We have studied the consequences of Bose-Einstein condensation of light ultrarelativistic particles when chemical freeze 
out occurs well before kinetic decoupling at Td < T c with T c the critical temperature below which a non-vanishing 
condensate fraction exists. We find that the presence of the condensate hastens the onset of the non-relativistic regime 
and that Bose-Einstein condensed particles can effectively act as a CDM component even when they decoupled being 
ultrarelativistic. The reason for this unusual behavior is that the particles in the condensate all have vanishing velocity 
dispersion. 

For thermal relics we find 
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The combination of data in ref. [14| from dSphs when applied to light thermal relics yields the mass range 
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with the implication that if these particles are suitable DM candidates, they must decouple at high temperature 
when the effective number of ultrarelativistic degrees of freedom is gd > 100. Namely, in absence of a BEC, thermal 
decoupling must occur above the electroweak scale. In the BEC case, for Td <C T c , the fulfillment of the bound 
requires very large gd- Namely, in the presence of a BEC thermal decoupling occurs at a scale much larger than the 
electroweak scale for <§; T c . 

Assuming that the DM particle is the only component with the density pdm today, we obtained an independent 
bound from velocity dispersion which for the favored cored profiles jl4j yield the lower mass bound 
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For light thermal relics this bound implies that to > 0.6 — 1.5 keV with a suppression factor Td/T c in the BEC case. 
For light thermal relics that decoupled while ultrarelativistic we find the primordial phase space density 

v , 3 / \ 3 ( 0.177 Fermions 
EEK „ 10 « eV / cm 3 gd ) 0.247 Bosons without BEC (6.6) 

° DM fkm/s) \ kcV / I 0.247 {T c /T d ) I Bosons with BEC . 



An enhancement factor (T c /Td)z appears in the r.h.s. in t he p resence of a BEC. 
For wimps with kinetic decoupling temperature lOMeV [33j |. we find 



°lim P A_/^ 3 I 100 GeV 



The observational data compiled in ref. [l4j assuming a favored cored profile suggests 
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If the distribution of dark matter is cusped, ref. [IH gives the value for the density p s ~ 2 TeV/cm 3 yielding 

^ ~ 2 x 10 9 eV/cm3 3 . (6.9) 
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Therefore, for gd > 10 the primordial phase space density for thermal relics with to ~ keV favors a cored distribu- 
tion. 

Notice that a bosonic thermal relic that features a BEC can behave as CDM with small velocity dispersion and a 
primordial phase space density consistent with cusped distributions if Td <C T c . However, these BEC DM candidates 
must decouple at a temperature scale higher than the electroweak. 

Recent results from ./V-body simulations suggests that the phase space density relaxes by a factor ~ 10 2 during 



gravitational clustering for < z < 10 [44|. Combining these numerical results with the observational results on dSphs 
[14| and the present DM density, we conclude that the mass of thermal relics that decoupled when ultrarelativistic is 

2 8 
m CO red ~ — keV , m cusp ~ — keV . (6.10) 

9 4 3 4 

The decoupling temperature for the DM candidate that would favor cusped profiles must be near a grand unified scale 
for a large symmetry group with gd > 2000 which effectively results in a colder relic today with a far smaller velocity 
dispersion. 

The enormous discrepancy between the primordial phase space density for WIMPs of to ~ 100 GeV; Td ~ 10 MeV, 
eq. (|6.7p and the phase space densities in dSphs, either cored feq |6.8l ) or cusped fea l6.9p cannot be explained by the 
two orders of magnitude of gravitational relaxation of phase space densities found with recent iV-body simulations 
p3 | , although these initialize the simulation with much smaller values of the primordial phase space density. 

We have studied a scenario for decoupling out of equilibrium motivated by previous studies of particle production and 
thermalization via an UV cascade. The distribution function obtained from previous studies [46j, remarkably describes 
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the non-equilibrium distribution functions for sterile neutrinos produced either resonantly fl9j or non-resonantly [l8j as 
well as a recently proposed model for halo structure [26| . Our bounds in terms of arbitrary distribution functions lead 
to the following bounds on the mass, phase space density and velocity dispersion of these light relics that decoupled 
out of LTE: 



• For sterile neutrinos produced non-resonantly via the Dodelson-Widrow mechanism [18j | we find 
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The upper and lower bound on the mass can only be compatible if the sterile neutrino decouples with gd > 20—30. 
For m ~ keV the primordial phase space density is compatible with cored but not with cusped profiles in the 
dShps data [l4| • Combining these bounds with the results from TV-body simulations on the relaxation of the 
phase space density [44| and with the observational constraint eq. l|4.30[) [14j |. we obtain the value 
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for the mass of sterile neutrinos produced non-resonantly by the Dodelson-Widrow mechanism. 
• For sterile neutrinos produced by a net-lepton number driven resonant conversion [lf| we find 
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The small velocity dispersion is a consequence of the distribution function being skewed towards small momen- 
tum. Again for to ~ keV, the primordial phase space density is compatible with cored but not cusped profiles in 
the dShps data [l4| . For sterile neutrinos produced by resonant conversion, a similar analysis as for the previous 
case yields 
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For the model proposed in ref. [26j we find 
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It is noteworthy that the TV-body results of ref. 44j which yield phase space relaxation by a factor ~ 10 2 bring 
the values of the primordial phase space density of the above cases within the range consistent with the phase space 
densities for cored profiles in dSphs [1J| for to ~ keV. On the contrary, in the case of WIMPs with to ~ 100 GeV, Td ~ 
lOMeV, relaxation by many orders of magnitude is necessary for their phase space densities to be compatible with 
the observed values both for cores and for cusps. 

Therefore the bounds eqs. (|6.12[) - (|6.14p confirm that ~ keV relics that decouple out of equilibrium while ultrarela- 
tivistic via the mechanisms described above yield values for phase space densities that are in agreement with cores in 
the DM distribution. 

The results obtained in this article for the new mass bounds, primordial phase space densities and velocity dispersion 
in term of arbitrary, but homogeneous and isotropic distribution functions establish a link between the microphysics 
of decoupling and observable quantities. They also warrant deeper scrutiny of the non-equilibrium aspects of sterile 
neutrinos [61| for a firmer assessment of their potential as DM candidates. 
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